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1. Introduction 


1.1 General 

The safety and reliability of structures has always 
been a matter of vital concern to the aerospace industry. 

In this respect, fracture mechanics (FM) is an especially 
useful technology, since it can provide a quantitative 
description of the capability of structural parts to 
tolerate flaws. Initially, FM concepts covered quasi-linear 
elastic conditions (LEFM) . Later, these methods were 
further developed to cover more general situations: a need 

existed to extend LEFM concepts to include cases where 
yielding was not necessarily contained in very small 
regions. This is often the case for new and tougher 
materials subjected to higher loads and used in thinner 
sections. This led to the development of the Elastic 
Plastic Fracture Mechanics (EPFM) Methodology. 

To apply these methods, two pieces of information are 
needed: the so-called material /specimen response to 

deformation, and the material response to crack extension. 
The former, obtained by finite element analysis or 
experimental calibration, consists of two expressions 
connecting the J-integral, load P, load-point displacement 
v, and crack length a for the specimen geometry of interest; 
the latter consists of a characterization of the way the 
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material resists crack extension for the type of load 
applied: J, or a similar parameter, versus crack extension, 

for monotonic load, da/dN versus AK or AJ for cyclic 
loading, or da/dt versus K, C* or Ct for creep crack growth. 
Within certain limitations, these curves are assumed to be 
specimen geometry independent, i.e., the curve obtained from 
a small laboratory specimen applies to the structural part 
under consideration, as well. 

A simple computer program can be developed to combine 
the two pieces of information mentioned and assess the 
reliability of the structural part of interest. 

It is very important to verify that the curve of 
material response to crack extension is geometry 
independent: the limitations of the parameters and/or 

approaches used must be understood, their validity limits 
must be clearly identified with the goals of improving the 
characterization of the phenomenon and proposing new 
parameters and methods to extend the range of applicability 
of existing models. 

1.2 Elastic Plastic Fracture Mechanics 

Specifically, for the case of EPFM applied to monotonic 
load, the mentioned limitations are expressed in terms of 
the amount of crack extension to ligament ratio, r, the 
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ratio of ligament to applied J over the yield strength, m, 
and the ratio of logarithmic increase of J to logarithmic 
decrease in ligament ,co. 

To overcome some of these limitations, particularly the 
one on r, Ernst [1.1] proposed a modified version of J 
called J M . Resistance (R) curves plotted in terms of J M 
were not subjected to the same limitations as those using J, 
and in general, showed a better correlation between 
specimens of different size and geometry. 

More recently, this methodology was further extended, 
general formulas were developed for Jm and Jp for growing 
cracks, criteria were proposed to identify the limits of 
applicability of both parameters, methods were presented to 
make use of the information of experimental points beyond 
this limit, and several schemes were proposed to extrapolate 
small laboratory specimen resistance curves to large amounts 
of crack extension, using Jy* Jd* or other parameters [ 1 . 2 ]. 

Although the progress made has been significant, and 
understanding has been gained on how to represent the R 
curve [1.2-1. 6], there are still several very important 
points that need to be addressed before the method can be 
safely applied. 

Among the most important ones, is the need to extend 
this whole methodology to include cases involving three 
dimensions (3D) : it is essential to know how specimen 
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thickness, constraint, and the possible dependence of the 
fracture mechanism on specimen thickness may affect the 
fracture resistance. 

Ultimately, this knowledge gained from "2D" planar 
specimens should be used to explain and predict the behavior 
of real 3D defects found in structures, e.g., surface or 
embedded cracks. 

1.3 The Leak Before Burst (LBB)_ Criterion 

Pressure vessels containing surface flaws are often 
required to comply with the so-called LBB criterion. LBB is 
understood as the condition where, upon application of a 
load, an existing flaw is assumed to grow through the wall 
of a pressure vessel in a prescribed fashion, causing 
leakage before the pressure vessel bursts. 

Pressure vessels of interest to NASA must comply to 
MIL-STD-1522A Standard General Requirements for Safe Design 
and Operation of Pressurized Missile and Space Systems. 

This document requires that: (1) the aspect ration, A.R. = 

a/2c (crack depth to total width ratio) needs to be in a 
range from 0.05 and 0.5, and (2) LBB will occur if 

*S/ > 2aB w withao < and a >1. (l.i) 

/ a cp * * 
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The rationale behind this expression is that the 
initial semi-elliptical flaw will grow in a self-similar 
manner, i.e., keeping A.R. = a/2c constant until the crack 
depth a is exactly equal to the thickness B, as shown in 
Figure 1.1. At that time, it is considered that the flaw 
becomes a through crack with a total length of 2c, with 2c - 
B/(A.R.) 0 , where the subscript 'o' is used to designate the 
initial state. Finally, to prevent the crack from running 
unstably in the longitudinal direction, i.e., bursting, the 
applied K given in equation 1.1 is required to be below the 
toughness Kj^ associated with this final flaw configuration. 

K = U- 2 ) 

*pp °p 

The MIL- Standard is inadequate for use with modern 
materials and designs because: (1) the above equation only 

holds for an aspect ratio A.R. = 0.5, (2) the flaw shape is 

considered to always remain elliptical with constant a/2c, 
and (3) the whole analysis is based on LEFM concepts. 

On the other hand, for the real materials, thicknesses, 
and typical flaws of interest, the situation is markedly 
different from the one assumed, above. A growing crack, 
shown schematically in figure 1.2, is marked with the 
numbers 1, 2, and 3 to suggest the evolving shape that is 
observed. Cracks do not grow in an elliptical self-similar 
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manner, but rather in a very complex shape, with a dimension 
in the direction parallel to the surface that is longer in 
the interior than on the surface. Moreover there is no 
guarantee that this dimension can be conservatively 
estimated by taking the original (2c/a) 0 and multiplying by 
B. 

1.4 This Pro-iect 

The EPFM Methodology has evolved significantly in the 
last several years. Nevertheless, some of these concepts 
need to be extended further before the whole methodology can 
be safely applied to structural parts. Specifically, there 
is a need to include the effect of constraint in the 
characterization of material resistance to crack growth and 
also to extend these methods to the case of 3D defects. 

As a consequence, this project was started as a 36 
month research program with the general objective of 
developing an elastic plastic fracture mechanics methodology 
to assess the structural reliability of pressure vessels and 
other parts of interest to NASA which may contain flaws. 

The project is divided into three tasks that deal with 
(1) constraint and thickness effects, (2) three- 
dimensional cracks, and (3) the Leak-Before-Burst (LBB) 
criterion . 
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Task 1 . Constraint and Thickness Effects 


This task includes the study of the problem of 
constraint and thickness effects observable in the 
fracture process, for different specimen sizes and 
geometries and for materials of interest. 

Specifically, the following subtasks will be 
performed: 

a) The large body of available data from centers around 
the World will be gathered to study this effect in 
specimens of different size and geometry. 

b) Resistance to crack growth tests will be conducted 
using specimens of different size and geometry, on at 
least one material of interest. The material will be 
provided by NASA; Georgia Tech will machine the lab 
specimens . 

c) Characterization of fracture surfaces to determine 
mechanisms of fracture, and typical surface 
dimensions will be performed using modern 
quantitative metallographic techniques. 
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d) Using the information obtained, models will be 

developed to describe the effect of constraint on the 
growth of cracks under elastic plastic conditions. 

Task 2. Three Dimensional Cracks 

The problem of applicability of EPFM concepts to 3D 
crack problems, in materials of interest, will be studied in 
this task. 

Specifically, the following subtasks will be performed: 

a) Plates containing surface cracks with different 
initial crack aspect ratios and relative crack-to- 
plate geometry will be tested. The evolution of the 
crack shape (planar) and the crack surface 
displacement with loading will be determined. 

b) Analytical and numerical efforts will be devoted to 
determine values of J and constraint along the crack 
front . 

c) The models and information obtained from Task 1 will 
be used here to predict the behavior of these 3D 
cracks . 

d) Predictions and experimental results will be compared 
and, if necessary, refinement of the models will be 
made. 
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Task 3. Leak Before Burst (LBB) Criterion 


The body of information obtained in the previous tasks 
will be organized in the spirit of the current LBB criterion 
into a Methodology format for the assessment of the 
structural integrity of parts containing defects. 

1 . 5 This Report 

This report covers the activities of the period March 
1994 through August 1994. In this period, full advantage 
was taken from the experience and knowledge gained in 
previous projects [1.6-1. 9]. In particular, some efforts 
were devoted in this project to complete and extend 
previously obtained results. 

In Chapter 2, experimental efforts to characterize 
three dimensional aspects of fracture present in "two 
dimensional", or planar configuration specimens have been 
continued. Chapter 2 specifically contains the discussion 
associated with the determination of, and use of, crack face 
separation data. 

In Chapter 3, the results of the fracture resistance 
testing for a variety of specimen configurations (J M R- curve 
format) is presented. The discussion presents the bases for 
like constraint for two materials of interest to NASA-MSFC 
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(in particular, 6061-T651 aluminum alloy and IN718 STA1 
nickel-base super alloy were characterized) in terms of the 
ligament dimensions, and compares these bases to the 
resulting Jj^R-curves. 
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^Assumed Crack Shape at Leak Condition 
/r- Assumed Self-Similar Crack Growth Steps 



Figure 1.1 Simple Leak Criterion 



Figure 1.2: Actual Growth of a Part-Through Crack 
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CHAPTER IX 


EEFQRIS IQ CHARACTERIZE THREE-DIMENSIONAL EFEfiCTS 
OBSERVABLE IN XttQ DIMENSIONAL FRACTURE SPECIMENS. 

by D. Lambert and H. Ernst 

2.1 Introduction 

The underlying purpose of this research is to develop a 
methodology which would allow the characterization of three- 
dimensional (3D) effects in fracture. This characterization 
should include: 

a) elastic plastic fracture behavior 

b) geometric effects arising from crack front 
curvatures (curvatures are present, for example, 
with surface cracks) , 

c) geometric effects related to thickness and 

ligament length (the gross sizing details that 
affect the three-dimensionality of the state of 
stress at the crack tip) , and 

d) loading geometry effects (including three- 

dimensionality of the far field stress arising 
from the character of applied loads, and, 



especially the gradient of the far-field stress 
arising from differing ratios of bending-to- 
tension) . 

Two-dimensional (2D) or planar specimens have been 
observed to generate different fracture resistance curves 
when different thicknesses are tested. Specifically 
discussed here are J^R curves that use the J-modified 
parameter as developed by Ernst [2.1 and 2.2] . Different 
configurations have been shown to support a differing degree 
of triaxiality of the stress field in the vicinity of the 
crack front, where the fracture process is occurring. The 
degree of stress field triaxiality that is exhibited is 
referred to as the constraint. The curves are a result 

of the different, averaged constraint in each specimen. 

Even though the configurations are considered to be planar, 
curvatures can develop in crack fronts that result from 
in the presence of a gradient of the constraint 
within the specimen. Thus, the complexities that occur in 
the most general cases of fracture appear in the simplest 
cases of planar specimens. Ultimately, to evaluate fracture 
resistance retaining a planar analogy requires that the 
crack front fall within specific limits of straightness. 

Since, 3D stress fields are present in planar 
configurations, an effort to map the crack face separation 
profiles of a variety of geometries as a function of the 
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position within the cross-section has been proposed. The 
crack tip opening displacement (CTOD) is a linear function 
of the J-integral [2.3]. Profiling represents an extension 
of that functional relationship. 

One goal of the overall research program is to test a 
wide variety of planar specimens, varying the thickness and 
length of the initial remaining ligament, as well as the 
bending-to-tension ratio due to the nature of the applied 
load to produce significant changes in fracture behavior. 
The results will be compared using the ligament dimensions 
as variables. Ultimately, the approach is expected to 
produce parameters and fracture behavior that can be 
generalized to the 3D cases that are of the most interest. 


The displacement of the faces of a crack are a function 
of the loading and of the position. Using a two dimensional 
analogy, the displacement of a point along a crack face 
within an elastic body was given by Tada, et al [2.4]: 



( 2 . 1 ) 
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Here, v e i is the elastic displacement in the loading 
direction at position r, measured from the crack tip to the 
point in question, is the equivalent modulus (E =E for 
plane stress, E'=E/(l-v 2 ) for plane strain, E is Young's 
Modulus, and v is Poisson’s ratio). The loading is 
specified in the presence of a flaw by K, the stress 
intensity parameter. This displacement relationship has a 
square-root of r form. 

Hutchinson [2.5] and Rice and Rosengren [2.6] developed 
a similar form for plastic response that follows Ramberg- 
Osgood deformation characteristics, i.e.: 

— = a(— r (2.2) 

a. 


In this equation, e and C are the equivalent strain and 
stress and e 0 , o Q and n are material constants. The form of 
the plastic displacement is, as follows: 


v p , = k-J$"r*" (2.3) 

This equation is written for a non-growing crack, and v p i is 
the displacement of the body, assuming Ramberg-Osgood type 
deformation, and k includes the functionality with regards 
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to the constraint, i.e., plane stress or plane strain 
condition. Looking at the equations (2.1) and (2.3), the 
constraint appears in the coefficient E' for the linear 
elastic case and in the coefficient k for the plastic case. 
Thus, the separation at various points through the ligament 
thickness could be expected to reflect that difference in 
constraint that arises with the depth into the thickness. 

It may also provide a measure of that constraint. 

Since the development above is for a non-growing crack 
situation, differences that occur between the theoretical 
elastic plus plastic displacements and the displacement 
profile of an actual growing crack near the crack tip might 
provide a fracture criterion on that local level. 

2.3 Profiling Matrix and Details 

One objective in the research was to characterize the 
degree of separation between the mating surfaces of cracks. 
This separation profile is a function of the level of J and 
of the position within the ligament. In this case the 
position would include the distance from the load-line, x, 
in the direction of crack growth, and the depth beneath the 
surface in the thickness direction, z. 

The primary effort in the past six months has been to 
analyze the crack face separation data of selected 
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specimens. The specimen identities and the corresponding 
configurations appear in Table 2.1, at the end of the 
chapter. Data was generated for a total of six compact 
tension (CT) specimens and three center-crack tension (CCT) 
specimens. The CCT configuration produces two crack 
profiles per specimen, and thus the total number of crack 
fronts observed is thirteen. Profiles were made of the 
AL6061-T651 only, because the IN718-STA1 material has proven 
too hard to polish in the same fashion as the aluminum, 
epoxy infused into the gap of the crack was effective in 
producing a well-defined crack profile for the aluminum, but 
the profiles of the nickel were rounded and poorly-defined. 
Until the techniques have been modified to overcome the 
rounding, the profiling was suspended for the nickel. 

After mounting the specimens in epoxy, the exposed 
surface was polished to provide a surface profile. After 
recording the profile, 0.025- to 0.035-inches was removed by 
grinding and polishing to produce the next profile to be 
recorded. This was continued for each specimen into the 
center of the cross-section. 

2.4 Profiling Results 

The crack tip opening displacement, 85 , was shown to be 
proportional to J or J M by Hellmann and Schwalbe [2.7] and 
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Ernst, et al [2.8], and is used as a fracture parameter. 
The crack grows away from the location where 85 is measured, 
and the measurement is expected to eventually lose 
sensitivity. It was thought that the shape of the loaded 
crack might provide a useful modification of the 85 
parameter, since the measurement location is being 
continuously updated as the crack grows. This section 
presents the results of efforts to establish the shape of 
the crack cavities of the compact tension and center-cracked 
tension specimens. The profiles are the "xy" plots of the 
crack shapes that were exposed by sectioning at different 
levels through the thickness. The results of the profiling 
will be given here. 

The information shown herein is organized as follows: 

( 1 ) observations of crack profiles through-the-thickness, 

( 2 ) comparisons of crack face separation profiles of 
selected specimens, (3) stretch analysis and the 
correlations, and (4) initiation stretch, and (5) a 
discussion of the observations. 

2.4.1 crack Profiles 

Although the individual profiles were not included, 
certain observations were made of their character, and these 
are summarized, below. 
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Figures 2 . 1 show the change of the crack shape as a 
function of the section depth for two specimens of different 
configurations. These figures show the shapes of the cracks 
at each of the section depths. Cross-sections were taken at 
each 0.025- to 0.030-inches starting at the specimen 
surfaces and continuing into their center planes. The 
fraction of depth below the surface is indicated by Z s 2z/B 
with Z = 0 at the surface and Z = 1 at the center plane. 
The vertical direction is the loading direction, and the 
horizontal direction is the direction of crack growth, with 
dimension referenced to the load line. The initial flaw 
location is indicated by "a 0 ". In figure 2.1a, the specimen 
is a CT specimen of thickness B = .85-inches (specimen #C9, 
CT, W = 2-inches, B = 0.85-inches, a/W = .5). At the 
surface, the crack profile is sharply angled beyond the 
fatigue precrack. This angled crack shape is associated 
with shear (note that the two axes of the figures are not of 
equal magnification, thus the angle shown is not the true 
angle) . For progressively deeper sections, the profiles 
become progressively straighter and longer. 

Figure 2.1b is a similar development to figure 2.1a. 
In this case, the specimen is identical except for the 
thickness of 0.25-inches (specimen #51, CT, W = 2-inches, B 
= 0.25-inches, a/W = 0.5). This set of crack profiles again 
shows the same angled character at the surface and the 
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tendency toward straightness as the center section is 
approached. 


2.4.2 Crack Face Displ acement Profiles 

An evaluation of the "calibration" value of the crack 
face displacement was computed as the sum of the elastic 
displacement plus an estimate of the plastic displacement. 
This total calibration displacement is a fitted function 
that agrees with the fatigue precrack region. The best fit 
to the precrack region appeared when the calibration was 
fitted using the length of the crack at the deepest point 
and had a form, as follows: 




c pl r; 


(2.4) 


In this equation, the calibration displacement is 6fi t (r), K 
is the average stress intensity factor, E* is the effective 
modulus, and r s h is the position referenced to the deepest 
point of the crack front in the whole cross-section. C p i 
and q are constants fitted to the precrack region. This 
calibration displacement was used to determine the stretch, 
discussed below. The stretch is calculated as the 
difference between the fracture profile and the calibration 
value. 
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Figures 2.2 - 2.4 show sets of the crack displacement 
profiles for five specimens. The precrack region has been 
omitted in the figures. Figure 2.2 presents one set of 
crack face displacement data for the entire cross-section of 
a typical specimen from surface to center. In this case, 
the specimen (#81, CT, W = 2-inches, B = 0.5-inches, a/W = 
0 . 75 ) is the same as the two earlier specimens, except for 
having a longer initial crack length and an intermediate 
thickness. The horizontal axis is the crack extension from 
the initial average crack tip location. The initial crack 
tip is at the zero location on the horizontal and at the 
beginning of the test, the displacement was zero on the 
vertical. As the load was applied, the initial crack tip 
blunts upwards and eventually, the crack begins to grow to 
the right. At zero crack extension (left side of the graph) 
is the displacement associated with the tip of the fatigue 
precrack at the final position of the test . The crack has 
grown to the right while the points to the left of the 
evolving crack tip have displaced upwards. The surface 
cross-section plots (lowest plots in the figure) show a 
substantial stretching or blunting before the crack begins 
£ o gjf 0 ^. This is seen as a vertical excursion (stretch) in 
the displacement profiles with no significant horizontal 
excursion (crack extension) . For example, the blunting of 
one of the surface plots is the segment labelled "AB" in 
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figure 2.2. For successively deeper sections, the stretch 
is progressively less, while in the cross-section the crack 
is progressively longer. At the center of the specimen (top 
plots in the figure), no discrete stretch zone is obvious. 
The data presented in figure 2.2 shows substantial noise, 
and so the plots have been filtered by a simple averaging 
technique to remove some of the noise without grossly 
affecting the shape of the displacement profiles. 

Each of the graphs in figures 2.2 - 2.4 show all plots 
associated with one of the specimens. Following the curves 
from the top-left to the bottom-right, once the initial 
stretching has occurred, all of the profiles have roughly 
the same shape, i.e., all of the displacement profiles have 
the same displacement at the location of the initial crack 
tip (crack extension of zero), and eventually, they develop 
approximately the same slope. The difference is that the 
surface profiles first stretched and then began to fracture, 
while those beneath the surface stretched less and the crack 
has grown longer. It should be noted that a few of the 
surface profiles may appear identical. This is also true of 
some of the center profiles. A set of 'transition profiles* 
can be seen that blends smoothly from the stretch— and 
fracture behavior of the surface to that of the center. 
Figures 2.3 show two identical specimens from a multi- 
specimen test subjected to two different input displacements 
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and exhibiting intermediate and long crack growth (CT, B = 
1/2-inch, a/W = 0.5, intermediate growth: #84, long growth: 

#E2 ) . For specimen #84, the surface displacements seem to 
be developing still, although the profiles are taking on the 
character of the early part of the profiles of specimen #E2 . 

The specific behaviors observed in figures 2.2 - 2.4 

were separated onto figures 2.5 - 2.7 to better illustrate 
certain observations. Figures 2.5 show the first three 
surface profiles for two specimens (#C9, CT, B = 0.85- 
inches, a/W = 0.5, and #D6, CCT, B = 1/2-inches, a/W = 0.5). 
In each graph, the data plotted together are quite similar. 
The Z = 2z/B = 0 plots are actually below the adjacent near- 
surface plots, but are quite similar. The crack has been 
grown to such a great length and substantial tunnelling has 
occurred, this is thought to be a secondary behavior arising 
from the thin surface section remaining once the center- 
plane crack has developed significantly. Figures 2.6 show 
that the center profiles for the three specimens stabilized 
to the same length and shape. Figures 2.7 show the 

transition profiles of one CT and one CCT specimen. In each 
case, if surface or center profiles showed the identical 
behavior, one of these was included for reference to the 
transition region, and thus the full range of profiles is 
represented. The profiles are taken (generally) at every 
.025 to .030 inches, and in figure 2.7a the profiles seem to 
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be equally spaced. This is not so obvious in the CCT 
specimen profiles in figure 2.7b. The author suspects that 
the behavior of the CCT profiles (figure 2.7b) would be best 
described as a surface effect that occurs for all but the 
inner-most sections. 

The similarities and transition characteristics 
described in the presentation above will be discussed in 
more specific detail below. 

2.4.3 Stretch Analysis 

The deviation of the separation profiles from the 
calibration profiles is a result of the growth of the crack 
from the initial position to the final position as the load- 
line displacement is being applied. For surface sections, 
the difference is quite pronounced, but for center sections, 
the deviation is insignificant. One point to be raised here 
is that the anticipated calibration profiles were expected 
to be parabolic, but it was seen that the fit was concave- 
upwards in many cases. The plastic component of the 
calibration function, the second term in equation 2.4 has an 
exponent q of less than unity. This exponent is shown in 
the EPRI handbook [2.9] as q = l/(n+l), which for this 
material is q « 0.048. Clearly, this is not the case. The 
fit using q = 0.048 does not fit the precrack region well. 
It was finally decided that the specimen configurations were 
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such that the single-term calibration was not accurate in 
the precrack region, and perhaps some beam deflection effect 
was being pronounced in the profiles. In support of this, 
two thoughts are offered: (a) the fit functions used do 
seem to fit the fatigue precrack region well with few 
exceptions, and (b) in the case of the center plots of the 
specimen #C9 (CT, W = 2-inches, B = 0.85-inches, a/W = 0.5), 
the concave upwards curve yields a similar shape to that of 
the data all the way to the tip of the crack. The use of 
q * 0.048 would require that the apparent crack length be 
shorter than the final crack. Further investigation into 
the proposed beam deflection effect or to find the 
appropriate calibration function might be warranted, however 
the results presented here are consistent and physically 
reasonable . 

The deviation between the calibration profile and the 
data is referred to as stretch. The next set of figures, 
2.8 - 2.12, show the stretch developed by the specimens at 
different section-depths. Figures 2.8 show surface stretch 
plots for four of the specimens. First, some of the stretch 
occurs devoid of crack growth, and then crack extension 
begins. Once the crack begins to extend at the section, it 
appears that the stretch may continue to increase, and this 
is true of all but the stretch profiles of the CCT specimen 
(figure 2.8d). The CCT surface stretch profiles, seem 
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horizontal. This increase may be a result of the use of 
improper calibration functions, although some inclination is 
expected. Otherwise, unstable fracture behavior would 
result. 

Moving to the center stretch profiles in figure 2.9, 
the initial stretch observed at the surface prior to crack 
extension is gone. The center stretch figure 2.9a and 2.9b 
are flat, suggesting that the stretch is not present, but 
only crack extension plus the presence of the calibration 
displacement . 

The transition from the surface stretch profiles to 
those of the center are shown in figure 2.10, these show the 
change of stretch characteristics as a function of the 
section depth for each specimen. The results presented in 
figure 2.10b should be taken very qualitatively, because in 
the middle sections of this recorded transition, the 
calibration profiles were found to fit best for the concave 
upwards calibration, thus these stretch profiles should not 
be compared, except at best to produce a trend. 

Figures 2.11 show all of the stretch profiles for two 
of the specimens. In figure 2.11a, the intermediate crack 
growth specimen (#84, CT, B = 1/2 -inch, a/W = 0.5) is 
presented. The figure shows a surface-to-center transition 
except for the two centermost sections, perhaps. No summary 
figure has been provided, but the behavior is identical to 
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that of the early growth portion of the long crack growth 
specimen (#E2, CT, B = 1/2-inch, a/W = 0.5). 

Figure 2.11b shows all of the stretch profiles for a 
thin specimen (specimen #51, CT, W = 2-inches, B = 0.25- 
inches, a/W = 0.5). Only the center-most two profiles 
appear to be similar. The observed stretch profile behavior 
is better characterized by a continuous transition from 
surface behavior to center behavior with no region of 
similar stretch characteristics. 

Figure 2.12 was included to show a comparison of the 
surface stretch characteristics of all of the specimens 
profiled. In this case, three behaviors appear, and these 
align by the different Rl ratios of the specimens profiled. 
The uppermost curves are of Rl = b 0 /B = 2, and both are 
from the same specimen. These curves show an increasing 
slope after the initial stretch phenomenon. The R L = 1 
surface sections produce a second, shallower slope, and six 
curves are present from two different specimens. The Rl = 4 
and the CCT surface stretch plots possess a surface stretch 
character that appears to be similar and is horizontal after 
the initiation stretch. For all of the specimens shown, the 
initiation stretch is approximately the same. The slope 
following initiation could be representative of different 
stretch-versus-crack extension characteristics for the 
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different configurations, but is more likely due to the 
choice of an inaccurate calibration function. 

No graphs of the center stretch profiles have been 
included, because it is thought that these would be 
misleading. The behavior can be summarized, as follows. 

a) The center stretch profiles show negligible initial 
stretch behavior before the crack begins to grow. 

b) Once the crack begins to grow in a given center 
cross-section, the stretch continues to evolve as the 
crack grows for the Rj_, = 2 and the CCT specimens 
(figures 2.9c and 2.9d), while for the Rl = 1 and the 
r l = 4 specimens it appears that the the stretch may 
increase only slightly as the crack extends (figures 
2.9a, 2.9b, and 2.12). Again, only a qualitative 
evaluation is recommended in the absence of verified, 
correct calibration functions. Some quantitative 
analysis will be performed in the next chapter, 
however the results should be considered to have a 
low level of confidence. 

2.5: Initiation Stretch 

When a crack extends into a structural ligament, the 
newly created surfaces are displaced from each other. The 
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displacement has been found to be related to the J integral 
[e.g. , 2.8]. 



Here, 8 is the crack tip-opening displacement, J is the 
value of the applied J-integral, Oy S is the yield strength, 
and d n is a parameter, approximately equal to unity, that 
accounts for the differences in constraint. The observation 
of the relationship between 8 and J has lead to the idea 
that a measured displacement might be used to describe 
fracture behavior. Various definitions for 8 can be 
offered, e.g., Hellmann, et al, [ 2 . 7 ] defined 85, the 
displacement of a 5-millimeter gage length positioned at the 
original crack tip. This measurement was taken between two 
gage points, although the crack extended to positions 
increasingly more remote from the measuring points. The 
profiling effort in this investigation extended the 
measurement methodology to the tip of the current crack 
position. A calibration displacement was approximated and 
the remainder of the total crack face displacement was 8 st , 
the stretch. This has been fit to a form as follows: 


(8 it -8 Itf ) = C It (x) p 


( 2 . 6 ) 
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The values of S re f and C s t come from a linear regression 
analysis of 8 st versus (X) P , with % being the crack 
extension (the distance in the profile from the initial 
crack tip position to the final crack tip position, X — ®f 
a Q ) . Whereas the value of C st is somewhat intuitive, that 
of 5 re f is not obvious: the regression analysis produced 

better correlation with 5 re f than if the regression was 
forced through 8 re f = 0. It will be noted that the exponent 
of % is P = 1/ (n+1) = .0483 « 0, and with values of C st « 0, 
it is suggested that the value of 8 re f might best describe 
the initiation stretch for a given profile. The data for 
the various fit parameters follow in the tables below along 
with a discussion of anomalies, special treatment received, 
and implications of the outcome. 

Several characteristics were observed in the profile 
fitting efforts that should be mentioned. First, although 
exceptions existed, the elastic and plastic displacement 
fits to the precrack region of the profiles produced 
relatively constant coefficients. For a given specimen, the 
coefficients were approximately the same for all of the 
profiles. One intended output to be extracted using the 
calibration displacement fit was the stretch. The resulting 
stretch traces were then fitted as a function of the local 
crack extension, % * af - a 0 , raised to a power of 

p = .0483. This corresponds to displacements associated 
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with power-law plasticity, however it should be noted that 
the stretch displacement is a deviation from the bulk 
material behavior, represented by the calibration 
displacement. The early part of the stretch curves seemed 
to be almost vertical, exhibiting blunting, and then the 
stretch became approximately constant. Thus, imposing the 
exponent of p seemed reasonable. The resulting coefficients 
from the linear fit of the stretch to X raised to p appear 
in each table, above. Two behaviors were apparent: (1) the 
values of ^-st were relatively constant for the surface 
profiles, and the center profiles were a different constant 
C st , and (2) the values of 6 re f seemed to be increasingly 
negative for the progressively deeper profiles. The form of 
the fit of the stretch in this material gives 

S.=5^ + C.(x) , -5«+CC.»6, (2. T) 

A factor of C has been inserted to produce the value of the 
second term in the middle expression in equation 2.7 at the 
knee of the curve. The approximation on the right in 
equation 2.7 is reasonable in the event that the exponent p 
is much less than unity (for AL6061-T651, p = .0483), 
producing relative insensitivity of 5 s t with respect to x 
above a nominal level of %• The value of the sum, 8i, 
appears in the right-most column in each table and is 
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monotonically increasing for most of the compact tension 
specimens, and is relatively constant for the center-crack 
tension specimens. Ultimately, two different thoughts arise 
from these observations: first, the ligament of the CCT 
specimens must be more uniformily deformed throughout the 
whole of the fracture process, only exhibiting a small 
initiation stretch. Second, at least for this material that 
approximates elastic-perf ectly plastic behavior, the 
initiation stretch, 8^, appears to be a function of s, the 
depth beneath the surface. A strict definition of s is 
necessary at this point. The first profile was taken at the 
root of the necked region formed on the side of the 
specimen. This was assigned a value of z = 0 inches. 
Because of the necking phenomenon, the z = 0 surface is 
somewhat below the original surface, and thus, s was defined 
as z plus the neck-root depth. Figure 2.13 shows the 
initiation stretch, versus s for all of the profiles. The 
lines connect adjacent profiles in each specimen. The 
initiation stretch decreases as s increases, but the 
correlation between specimens is not as obvious. 

Various forms were considered to establish the nature 
of the behavior of 8j[ versus s, and clearly the best 
correlation was with a parameter sB/b 0 s k, and this 
correlation is plotted in figure 2.14. Several features 
visible in the figure must be discussed. 8^ increases from 
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the surface profile to some peak value and then decreases as 
k increases further. With the exception of the surface 
behavior and the behavior of the data indicated with dashed 
lines, all of the behavior seems to be along the same curve 
within experimental error. No obvious explanation exists 
for the results of the CCT specimen in the figure. Figures 
2.15 have been included to segregate certain behaviors 
observed. The surface profiles, the Rl = 4 data, and the 
CCT data points have been excluded from the data set in 
figure 2.15a. The initiation stretch data that remains in 
this figure shows a consistent trend. Figure 2.15b shows 
the surface points and includes all of the data from the two 
specimens with R L = 4 . The horizontal axis position is the 
location of the root of the neck expressed as K. The curve 
of the surface plots appears to be almost vertical, except 
that the lowest data points do show an increase in k from 
the other data. A possible explanation for the behaviors in 
the figures 2.15 is this: for each specimen, a surface zone 
exists, and as the depth below the surface increases, the 
surface effect makes a transformation into the uniform 
subsurface behavior that is shown in figure 2.15a. In the 
case of the two specimens with Rl = 4, the surface effect 
permeates the whole thickness, and the entire cross-section 
degenerates into a surface effect defined simply by the 
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calibration profile up to the point of initiation of crack 
growth. 

The maxima of each 8 i~versus-s curve occurs on the 
"uniform subsurface" curve shown in figure 2.15a, and the 
maximum initiation stretch can be estimated by the median 
value of 5iB/b 0 - .03 * constant. The minimum value is 

zero . 

The results suggest that given a method to evaluate the 
constraint along an arbitrary crack front, and the boundary 
conditions, the calibration displacement of each point could 
be determined, and an estimate of the stretch and crack 
growth characteristics of the point might be easily 
determined by comparing to the results of such a planar 
crack growth analysis as presented, above. 

2.6 Conclusio ns and Recommendations 

1. The profiling procedures used to generate the crack 
face separation throughout planar specimen cross- 
sections showed that the 85 -displacement, measured 
at the surface was the same as the displacements 
observed inside the cross-section. 

2. The 85 -displacement that is a common throughout the 
specimen cross-section at the initial crack tip is 
developed in a different manner depending on the 
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depth into the thickness considered The onset of 
crack extension was termed the initiation stretch, 
8^. The value of 5i seemed to be a function of the 
depth beneath the surface, although an alternate 
functionality that included the ligament 
proportions and the depth beneath the surface was 
also suggested. 

a) At some depth, the stretch prior to initiation is 
negligible. This is thought to be the plane 
strain condition. 

b) For all of the cross-sections, the crack angle was 
approximately constant through the cross-section 
after crack growth initiates. 

3. The use of profiling techniques shows some promise 
in the characterization, and profiling or some 
nondestructive technique (NDT) , such as computer 
tomography, might be employed for a different 
material, e.g., IN718-STA1, to compare with results 
of the AL6061-T651. The destructive nature of the 
profiling procedures is not expected to affect the 
results much, but NDT techniques would avoid the 
potential influences altogether. 
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2.8 Tables fo r Chapter 2 


Table 2.1: Matrix of Specimens Profiled 


(W, B , and b are in inches, a/W is nondimens ionall 


Spec 

Conf ig 

W 

B 

H 

a/W 

Remarks 

No of 
Profs 

E2 

1T-CT 

2 

1/2 


.5 

Baseline Specimen 

8 

81 

1T-CT 

2 

1/2 

.5 

.75 

Larger Init . a/W 

8 

51 

1T-CT 

2 

1/4 

1 

.5 1 

Thinner Section 

5 

C9 


2 

.85 

1 

.5 

Thicker Section 

13 

82 

1T-CT 

2 

1/2 

1 

.5 

Multi-Specimen 

8 

84 

1T-CT 

2 

1/2 

1 

.5 

Multi -Specimen 

8 

08 

1T-CT 

20%SG 

2 

1/2 

m 

.6 

Side-Grooved 

8 

D6 

CCT 

1 

1/2 

.5 

.5 

Tension 

2x8 

55 

CCT 

1 

1/8 

.5 

.5 

Tension, Thin 

2x3 

B9 

CCT 

1 

.85 

.5 

.5 

Tension, Thick 

2x9 


Table 2.2: Fitting Constants for CT Profiles 




Prof i 
le 

s 

Ar sh 


Cnl 

q 

8 ref 


8i 

Comments 

8201 

.01 

2 

-.216 

H 

.003 

.99 

X 

X 

X 



.03 

9 

-.204 

■a 

.002 

-.05 

X 

X 

X 

m **0 

8203 

G9 

-.192 

■■ 

.006 

.87 

X 

X 

X 


8204 

.09 

3 

-.175 

■i 

.007 

.95 

X 

X 

X 


8205 

.12 

5 

-.137 

■■ 

.009 


X 

X 

X 


8206 


-.080 

■■ 

.011 

.96 

X 

X 

X 


8207 

.20 

7 

-.053 

■a 

.006 

.86 

X 

X 

X 


8208 

.23 

8 

0 

■a 

.008 

Bgj 

X 

X 

X 
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Table 2.3: Fitting Co nstants for CT Profiles 


Prof . 

s 


" • 

f *0 

C-nl 

q 

mm 

warn 

*i - 

Comments 

8401 

.022 

-.237 

Bl 

.042 

.58 

-.14 

.20 

.039 

Cfit* 

large 

8402 

HSI 

D^i 


.035 

.71 

-.07 

.10 

.026 


8403 

.075 


1 

.034 

.70 

-.06 

.09 

.022 


8404 

■Uiia 

-.163 

1 

.033 

.67 

-.09 

.12 

.021 


8405 

.136 

-.118 

1 

.035 

.61 

-.13 

.16 

.016 



HUH 


1 

.033 

.65 

-.11 

.13 

.010 


8407 

■ren 

-.038 

1 

.034 

.64 

-.13 

.15 

.007 


8408 

.210 

0 

Bl 

.030 

.62 

-.02 

.02 

.004 

small 


Table 2.4: Fitting Co nstants for CT Profiles 


1 Specimen 

E2: AL, 

1T-CT, B=l/2 inch 

a/W= . 5 , Long 

Crack Gr< 

DWth 

Prof i 
le 

s 

Ar sh 

i 


q 

5 ref 


h 

Comments 

E201 

.020 

-.160 

H 

■Bl 

.83 

-.16 

• 21 

.029 

■KiitwHi 

E202 

.045 

-.325 


HHsBB 

.95 

-.24 

.30 

.031 


E203 

■tiiii 

-.308 

1 

mssm 

1.0 

-.30 

.36 

.024 


E204 

.101 

-.272 

1 

.114 

1.2 

-.34 

.38 

.012 


E205 

.134 

-.237 

1 

.114 

1.1 

-.26 

.30 

.013 


E2 06 

.161 

-.114 

1 

.114 

1.1 

-.36 

.40 

.002 


E207 

.192 

-.037 

II 

.112 

HI 

-.24 

• 26 

.001 

small 

E208 

.223 

0 

H 

.110 

HI 

-.16 

.18 

.002 

Km 
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2.9 Figures f nr Chapter 2 
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Figure 2.1: Typical cross-sections of cracked specimens, showing crack profiles from 

surface to center plane: (a) specimen #C9 (CT, B = 0.85-inches, a/W = 
0.5), and (b) specimen #51 (CT, B = 1/4-inches, a/W = 0.5). 
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Figure 2.2: Crack face separation profiles for cross-sections from surface to center 

plane of specimen #81 (CT, B = 1/2-inch, a/W = 0.75). Note: the 
normalized depth is 2z/B = 0 at surface and 2z/B « 1 at center. 
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Crack Extension, In. 



Figure 2.3: Crack face separation profiles for cross-sections from surface to center 

plane of two specimens from multiple-specimen test (CT, B = 1/2-inch, 
a/W = 0.5): (a) intermediate growth, specimen #84, and (b) long growth, 
specimen #E2. Note: the normalized depth is 2z/B = 0 at surface and 2z/B 
* 1 at center. 
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Figure 2.4: Crack face separation profiles for cross-sections from surface to center 

plane of two specimens of different configuration (a) specimen #51 (CT, 
B = 1/4-inch, a/W = 0.5), and (b) specimen #D6 (CCT, B = 1/2-inch, 
a/W = 0.5). Note: the normalized depth is 2z/B = 0 at surface and 2z/B = 
1 at center. 
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Figure 2-5: Crack face separation profiles for near-surface cross-sections of two 

specimens: (a) specimen #C9 (CT, B = 0.85-inch, a/W = 0.5), and (b) 
specimen #D6 (CCT, B = 1/2-inch, a/W = 0.5). Note: the normalized 
depth is 2z/B = 0 at surface and 2z/B * 1 at center. 
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Figure 2.6; Crack face separation profiles for center plane cross-sections of three 
specimens: (a) specimen #81 (CT, B = 1/2-inch, a/W = 0.75), (b) specimen 
#C9 (CT, B = 0.85-inch, a/W = 0.5), and (c) specimen #51 (CT, B = 1/4- 
inch, a/W = 0.5). Note: the normalized depth is 2z/B = 0 at surface and 
2z/B « 1 at center. 
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Figure 2.7 • Crack face separation profiles of two specimens from surface to center 
plane: (a) specimen #C9 (CT, B = 0.85-inch, a/W = 0.5), and (b) specimen 
#D6 (CCT, B = 1/2-inch, a/W = 0.5). Note: the normalized depth is 2z/B 
= 0 at surface and 2z/B « 1 at center. 
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Figure 2.8: Stretch-versus-crack extension data for near-surface cross-sections of four 

specimens: (a) specimen #81 (CT, B = 1/2-inch, a/W = 0.75), (b) specimen 
#C9 (CT, B = 0.85-inch, a/W = 0.5), (c) specimen #E2 (CT, B = 1/2-inch, 
a/W = 0.5), and (d) specimen #D6 (CCT, B = 1/2-inch, a/W = 0.5). Note: 
the normalized depth is 2z/B = 0 at surface and 2z/B « 1 at center. 
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Figure 2.9: Stretch-versus-crack extension data for center plane cross-sections of four 

specimens: (a) specimen #81 (Cl , B = 1/2-inch, a/W — 0.75), (b) specimen 
#C9 (CT, B = 0.85-inch, a/W = 0.5), (c) specimen #E2 (CT, B = 1/2-inch, 
a/W = 0.5), and (d) specimen #D6 (CCT, B = 1/2-inch, a/W = 0.5). Note: 
the nor malize d depth is 2z/B = 0 at surface and 2z/B *» 1 at center. 
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Figure 2.10: Stretch-versus-crack extension data for cross-sections from surface to 
center plane of four specimens: (a) specimen #81 (CT, B = 1/2-inch, a/W 
= 0.75), (b) specimen #C9 (CT, B = 0.85-inch, a/W = 0.5), (c) specimen 
#E2 (CT, B = 1/2-inch, a/W = 0.5), and (d) specimen #D6 (CCT, B = 1/2- 
inch, a/W = 0.5). Note: the normalized depth is 2z/B = 0 at surface and 
2z/B « 1 at center. 
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Figure 2.11: Stretch-versus-crack extension data for cross-sections from surface to 
center plane of two specimens: (a) specimen #84 (CT, B = 1/2-inch, a/W = 
0.5), and (b) specimen #51 (CT, B = 1/4-inch, a/W = 0.5). Note: the 
normalized depth is 2z/B = 0 at surface and 2z/B = 1 at center. 
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Figure 2.12: Comparison of stretch-versus-crack extension data for near surface cross- 
sections of various specimens: R L = 1, 2, 4, and CCT with R L = 1. 
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Figure 2.13: Comparison of initiation stretch as a function s, the depth below the 
original specimen surface for a variety of specimens. Note: at surface, s = 

0, and at center, s = B/2. 
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Figure 2.15: Comparison of initiation stretch as a function sB/b, the modified depth 
below the original specimen surface for a variety of specimens: 

(a) subsurface data (excludes first surface point, CCT data and R L = 4 
data), and (b) "surface data". 
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CHAPTER HI 


RESISTANCE CURVE ANALYSIS 

by D. Lambert and H. Ernst 
3.1 J and Resistance Curve Analysis 

In the literature, variations in the resistance curves 
of a given material have been reported for differently sized 
and shaped ligaments. This is due to the variations of 
constraint that occur for each configuration. To 
demonstrate the effects of constraint and for the purposes 
of this investigation, a test matrix was established that 
provided a wide variety of specimen configurations. The 
ranges of the configuration defining variables investigated 
appear in brief in Table 3.1. 

The different configurations tested produced different 
groups of resistance curves for both materials investigated, 
and the purpose of this section is to present these 
qualified results. A close look at the valid portion of 
each of the resistance curves would show two features: (1) 
the qualified J-modified resistance curve extends much 
further in most cases than that of the J-integral resistance 


curve, and (2) the J- integral and J-modified agree well up 
to the limit of J. It should also be noted that J-modified 
appears to be a smooth extension of the valid portion of the 
J-integral curve in each case. 

3.1.1 JR and J ^ R Curves for Compact Tension Specimens 

Figure 3 . 1 suggested that certain specimens should be 
compared, and the basis for comparison of the aluminum 
results was using Rl e b 0 /B, while B was found to be 
appropriate for the IN718-STA1. The next set of figures 
(figures 3.2 and 3.3) correlates the valid J and the 
qualified resistance curve data for the aluminum. The 

correlation has improved dramatically for both the J and the 
Jj 4 curves. 

Focusing on the Jj,j results at this point (figure 3.3) , 
the curves chosen for comparison are very similar. Notice 
that the Rl = 1/2 and the side-grooved configurations are 
shown in figure 3.3a, while those for Rl = 1 appear in 
figure 3.3b, and those for R l = 2 and 4 appear in figure 
3.3c. In figure 3.3c, a slight variation exists between 
curves for R l = 2 and R L = 4. The curve variation could be 
due to configuration (Rl ratio) differences. This will be 
discussed, below. Figure 3.4 has been included to show the 
improved correlations and the presence of different 
behaviors of the J M curves after qualifying the data. Three 
distinct curves of different slopes clearly are present in 
this graph. One of the side-grooved CT specimens shows data 
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1 curves (specimen #E3 , 


more reminiscent of the Rl = 
aluminum, CT, W = 2-inches, B = 1/2-inch, 20% side-grooved, 
a/W = 0.5): the side-groove was not directly on the 

precrack tip and the crack grew into the side-groove instead 
of starting there. The resulting crack front had 

considerable curvature, and a substantially higher 
resistance curve. Normally the side-grooved configuration 
produces a straight crack front, and this data is not 
thought to clearly fit into any particular category being 

investigated and will be ignored. 

Next, the "valid" J M data for the nickel alloy is 
presented in figures 3.5. A graph for the nickel alloy of 
thickness B = 0.80-inches was excluded at this point, 
because only one specimen was represented, and no comparison 
was available. The data segregated onto each figure 
compares well. Figure 3.6 compares all of the qualified J- 
modified data for the nickel alloy including the specimen 
omitted earlier (specimen #A5, nickel, CT, W = 2-inches, B = 
0.80-inches, a/W = 0.5) , and after the qualifying step, 

four distinct curves are apparent. The imposing of the Jm r 
inflection point limitation as the final point of 
consideration of the data has provided good correlation with 
the qualified data. 

The subsection that follows is devoted to a discussion 
of the results for the other two configurations used to 
analyze the fracture resistance of the AL6061-T651, namely 
center-crack tension and single-edge notch tension 
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specimens. This is followed by a section discussing the 
nickel CCT fracture behavior, and then a general discussion 
of all of the results. 

3.1.2 Aluminum rcT and s ent Fracture Res i stance Curves 

The AL6061-T651 center- crack tension (CCT) specimen 
data were shifted in the same manner as the CT data, 
validity limits were imposed, and the results appear in 
figure 3.7. The early part of the CCT specimen Jm r curves 
appear to follow the apparent blunting line to a much higher 
value of Jm than was observed with the CT specimens. It 
must be noted that the center-cracked tension specimen with 
side-grooves (specimen #C0, aluminum, CCT, 2W = 2 -inches, B 
= 0.85-inches, 20% side-grooved, a/W = 0.6) underwent 

unstable fracture, and the data was developed by assuming 
that the crack curvature correction factor remained constant 
throughout the test. The resulting resistance curve fits 
the trend of the other data in figure 3.7, and this is 
considered to be a reasonable estimate, since the side- 
grooved configuration is thought to produce a consistently 
straight crack front. The thinnest specimen (specimen #55, 
aluminum, CCT, 2W = 2-inches, B = 1/8-inch, shown as filled 
squares in figure 3.7) follows that of the thickest and the 
side-grooved specimens. This is not the behavior expected, 
but the data is thought to be bonafide. 

The aluminum single-edge notch tension (SENT) specimen 
test results are shown in figures 3.8. The first graph, 
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figure 3.8a shows the SENT data separately, while the 
second, figure 3.8b shows a comparison of the SENT data 
superimposed on that of the CT specimen data. The SENT 
results are the filled symbols in figure 3.8b and are very 
similar to those of the aluminum CT specimens. The only 
difference apparent is that the SENT specimen with ligament 
proportion R L = 1 follows the blunting line to a higher 
value of J than the CT specimen of the same Rl ratio. A 
closer scrutiny comparing the behavior of the aluminum for 
all three specimen types will be presented below. 

Figure 3 . 9 shows the data for all of the aluminum 
specimens on one graph. The blunting line has been shown 
for reference. In this figure, the CT and the SENT 
specimens of equal Rl value would not be discernable from 
each other within a small range of variation of the graphs, 
and the CCT specimens do not deviate from the blunting line 
until much higher values of J^. 

3.1.3: Nickel CCT Fracture Resistance Curves 

The IN718-STA1 center-crack tension specimen data are 
shown in figure 3.10. The behaviors of both specimens are 
very similar in the upper end of the J^R curves. The early 
behavior of the nonside— grooved specimen (specimen tt24, 
nickel, CCT, 2W = 2-inches, B = 1/4-inch, a/W = 0.5) has 
been attributed to the imprecision with which the crack 
lengths were measured. Figure 3.11 compares the various 
compact tension specimen results with those of the center- 
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Also shown in this figure is the 


crack tension specimens. 

ASTM blunting line for this material. As was seen with the 
results from the aluminum specimens, the nickel CCT 

specimens follow the blunting line to much higher J levels. 
In this case, the data does not appear to deviate from the 
blunting line for the duration of both tests. A more in 
depth discussion of the J M resistance curve results follows. 

3.1.4: Discus s ion of Jjj Resistance Curve Results 

Several characteristics that are very prominent in the 
adjusted curves should be discussed. The first comments are 
general, relating to both of the materials and focusing 
especially on the level of confidence in the data, on the 
success of J or in correlating the data, and on the 

behavior of the specimens. These comments are followed by 
characterizations of the fracture behavior for the two 
materials and a discussion of the differences in the 

perceived behavior between the two materials. 

The test data was rough, however this is not expected 

to invalidate the results or to reduce the impact of the 

conclusions drawn, herein. This roughness is thought to 
arise primarily due to imprecision associated with the 
determination of the intermediate crack lengths early in the 
test from measured unloading compliance values. The data 
still show definite trends, and all of the data is 
consistent with these trends with few exceptions. These 
exceptions will be mentioned at the appropriate time. 
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From the CT data, J^ seemed to produce results that 
were identical to the J-integral up to the established ASTM 
E1152-87 [3.1] limits of validity of J. A close look at the 
various JR- and J^R-curves would suggest that the point of 
deviation of the JR and J^R curves is very close to the same 
point indicated by the ASTM limit. One other possibility 
for the marker of the divergence point of J and Jj* was the 
maximum load. This occurred consistently at a crack 
extension of five-percent of the initial ligament length. 
The difference in J and J^ grows to about three-percent of J 
at the ASTM limit (i.e., Aa/b 0 = 0.1). By limiting J M data 
to include only data prior to the upwards inflection point 
visible in some J^R curves, the correlation of the data was 
improved throughout the data set with the only apparent 
variations being due to experimental error, and the J^R 
curve was a smooth extension of the JR curve beyond the 
limit of the J-integral. The physical phenomenon that 
becomes significant at the inflection point of J^ is 
believed to be the onset of large-scale plasticity, i.e., 
the saturation of plasticity into the back face of the 
specimen [3.2]. This subject was discussed earlier in this 
chapter. 

Using J^ resistance curve data, the fracture 
performance of the AL6061-T651 specimens is summarized in 
figures 3.3, 3.7, 3.8a, and 3.9. For the CT specimens, no 
specific thickness or ligament length effect is obvious in 
the data. Instead, the fracture resistance was a function 
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of the ligament length-to-thickness ratio, b 0 /B. Any 
secondary dependence of the fracture properties due to 
length or thickness of the ligament was overwhelmed by 
experimental error. Two sizes of the CT specimen 
configuration were used, and these were the 1/2T and the IT 
configurations, associated with planar dimensions defined by 
W = 1-inch and W = 2-inches, respectively. No difference in 
the results was found between these two configurations. The 
CCT and SENT specimens used the same initial ligament 
length, thus the ligament length, thickness, and proportion 
effects are not separable in the results. If expressed as a 
function of the ligament proportion, the trend of all data 
was identical to that of the CT specimens, i.e., the 
fracture behavior of the aluminum CT, CCT, and SENT 
specimens appear to be controlled by the ligament length-to- 
thickness (R l s b 0 /B) ratio. The observed behavior of the B 
= 1/8 inch CCT specimen (#55, aluminum, CCT, 2W = 2 -inches, 
B = 1/ 8-inches, a/W = .5) will be noted here as the one 
exception. 

Looking at figure 3.3a, aluminum CT specimens with R L = 
1/2 yields behavior identical to that of the plane strain CT 
specimens (side-grooved configuration). Next, looking at 
figure 3.3c, CT specimens with Rl = 2 and Rl = 4 appear to 
produce equivalent results to each other. Figure 3.4 shows 
all of the aluminum CT specimen results on one graph. The 
curves become more steeply inclined as Rl increases. 
However, the behavior of the J^R-curves seems to start from 
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a "saturated" value of plane strain that is present for the 
side-grooved and the Rl = 1/2 specimens, and then the 
curves appear to "saturate" again to a value that is taken 
to be plane stress when Rl £ 2. Thus, an active region of 
changing constraint for the aluminum alloy appears to be 
plane strain = 1/2 £ Rl £ 2 = plane stress. Within these 
limits, the average constraint is between plane strain and 
plane stress. 

The fracture performance of the IN718-STA1 CT specimens 
is summarized in figures 3.5 and 3.6 and figures 3.10 and 
3.11. The fracture behavior of the nickel alloy is not 
defined by the ligament proportion as is the aluminum, but 
is characterized by the ligament thickness, B. This is true 
for the CT specimens tested, but the nickel CCT specimens 
were limited to one thickness (B = 1/4 inch and B = 1/4 
inch, with side-grooves) , and no deviation from blunting 
line behavior was obvious in either case (figure 3.11). 

In figure 3.5a, the J M R curve for the nickel CT 
specimen with the longest ligament (specimen #14, nickel, 
CT, W = 2-inches, B = 1/4-inch, a/W = 0.5, i.e., b Q = 1) 
appears to be slightly higher that those for specimens with 
shorter ligaments, but this is thought to be associated with 
experimental error and can be easily influenced by the 
shifting discussed early in this section. Figure 3.6 
combines all of the nickel CT data to compare the results. 

The nickel J M R curves are found to incline more steeply 
as the specimens become progressively thinner. Figure 3.12 
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has been included to compare the results of the plane strain 
condition (side-grooved configuration) to those of the 
thickest specimen (specimen #A5, nickel, CT, W = 2-inches, B 
= 0.8 inch, a/W = 0.5). The axes are greatly expanded in 
this case, and the behavior of the B = 0.8 inch specimen 
appears to continue rising while that of the side-grooved 
specimens appears to be constant. A thickness of B = 0.80- 
inches should not be considered to be exhibiting plane 
strain fracture. 

The results of all of the CCT and SENT data are 
summarized in figures 3.7 - 3.9 (aluminum), and figures 3.10 
and 3.11 (nickel). The data for the CCT specimens was 
somewhat truncated, because unstable cleavage seemed to 
occur with only a small amount of crack growth. It should 
be noted that in all cases, the thickness was varied and the 
ligament length was kept the same at b D = 1/2 inch. 

The aluminum CCT specimen data appears in figure 3.7. 
With the exception of the thinnest specimen (filled square 
symbols, specimen #55, aluminum, CCT, 2W = 2-inches, B = 
1/8-inches, a/W = 0.5), the J M R cuves of the remaining 
specimens show different slopes, behavior congruent to the 
CT behavior: the specimen with side-grooves and the one 
with the lowest value of ligament proportion (Rl = 1/2) 
produce the lowest curves, while for increasingly higher 
values of Rl, the JR and J^R curves are increasingly higher. 
The behavior of the thinnest specimen (#55) was not as 
anticipated, however all of the other specimens provide 
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consistent behavior. It appears that the thickest specimen 
(specimen #B9, CCT, 2W = 2-inches, B = 0.85-inches, a/W = 
0.5) was sufficiently thick to produce plane strain 
behavior, however with the behavior of the thin specimen and 
with the sparsity of the data, conclusions cannot be stated 
as strongly as with the compact tension specimens. The 
aluminum CCT data does follow the same trend as that of the 
aluminum CT specimens, though, and the correlation of 
behavior with respect to R L appears similar. One large 
difference between the CCT and CT results is that curves for 
the CCT specimens diverge from the blunting line at a much 
higher value of than do the curves for the CT specimens 
(figure 3.8). 

The SENT specimen configurations were quite similar to 
some of the CT specimens. The SENT specimen has a large 
distance between the loading points, and the specimen halves 
will not rotate as the test proceeds to large displacements 
and long crack lengths. The results of the SENT specimens 
are similar to those of the CT specimens with similar R L - 
value, except that the SENT configuration continues to a 
higher J M - value before diverging from the blunting line 
(figure 3.7b). 

The nickel CCT specimens follow the blunting line to a 
much higher level of J M than the CT specimens before showing 
any deviation from the blunting line. In fact, the two 
nickel CCT specimens tested did not show any deviation from 
blunting line behavior for the limits imposed on the test. 
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The aluminum and the nickel show dramatic differences 
in fracture behavior. The aluminum specimens showed 
fracture behavior controlled by Rl = b 0 /B, and it appears 
that at Rl = 1/2, the behavior is equivalent to that of 
plane strain (side-grooved) specimens. If Rl is increased 
from 1/2, the JR and J^R curves become higher, up to a 
point. The curves of Rl = 2 are roughly equivalent to that 
of Rl = 4. This suggests that at Rl = 2, the ligament is 
subjected to plane stress fracture or the surface effect has 
saturated the ligament, and that fracture behavior will no 
longer change as Rl continues to increase. The nickel 
specimens were effected by thickness B, alone. 

The behavior of the nickel is quite different: over 
the range of compact tension specimen configurations tested, 
the fracture behavior of IN718-STA1 changes with thickness, 
only. Figure 3.13 suggests that the resistance curves for 
the IN718-STA1 compact tension specimens of ligament 
proportion Rl = 1 differ greatly. 

3 . 2 Power-Law Fit of R-Curve Data and Estimates 

For each specimen, the qualified JR and J^R data showed 
good agreement, although the limiting of the JR-curve to a 
crack extension of ten-percent of the initial remaining 
ligament length allowed a slight divergence between J and 
J M . The resulting R-curves showed certain common behaviors, 
and the aluminum and the nickel alloys exhibited different 
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characteristics: specifically, the aluminum data were 
identical for equal ligament initial length-to-thickness 
ratio, Rl s b 0 /B, while the nickel data aligned for like 
ligament thickness, B. The R-curves for the aluminum CT 
specimens exhibited an increasing slope as Rl increased, 
while the R-curves for the nickel inclined more steeply as 
the ligament thickness decreased. 

Little redundancy was included in the matrix of CCT 
test specimens and all of these specimens had the same 
initial remaining ligament length, thus thickness and Rl - 
proportion effects could not be separated. The aluminum CCT 
specimen data did produce increasing slope of the R-curves 
for increasing Rl (or alternately for decreasing thickness, 
B) , this is the same characteristic R-curve slope trend as 
the aluminum CT specimens. The nickel CCT specimens seemed 
to follow the blunting line, although the results of a 
regression analysis of the data has been included. 

The correlated data were fitted to a power law form: 

J M = J 0 -(Aar f 3 - 1 * 

In equation 3.1, Aa is the crack extension and J Q and m are 
fitted constants from a linear regression analysis of the 
relationship of the logarithms of Jjj and Aa. The fit was 
performed only for Jm because the common bases for 
comparison (Rl for aluminum and B for nickel) arose in the 
J M analysis, not in the J-integral analysis. The fitting 
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procedure used all of the qualified and correlated data in 
each case. The results of the regression analysis appear in 
tables 3.2 and 3.3, below. 

The data in table 3.2 generally show a consistent trend 
for both J Q and for m within the specimen types, except for 
the CCT specimen with Rj_, = 2 . The data for this specimen 
was observed to follow the blunting line, and the value of m 
= 1.01 appears to corroborate this. Generally, the values 
for the exponent m for the CT or for the CCT specimens 
should be considered to be relatively equal, except that the 
side-grooved configurations (marked "SG") have exponents 
that are noticeably lower than the nonside-grooved 
specimens. The crack extension will be large for a small 
increment of J M relative to the other configurations once 
the loading increases beyond the initiation value of Jm* 
The exponents, m, for the CCT specimens appear to be 
slightly higher than those of the compact tension specimens, 
but this may not be an appropriate conclusion: the 
difference may be a result of the fewer numbers of data 
points for the CCT, the relatively limited crack growth in 
the CCT specimen configurations, or some problem in 
providing an appropriate estimate of the average crack 
length. 

The regression analysis results for the nickel are 
self-consistent and follow the same trend as those described 
for the aluminum. The values of J D are much larger in this 
case, reflecting the higher toughness of the nickel. The 
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exponent m also appears to decrease more over the range of 
configurations than the aluminum. For B = 1/4-inch nickel 
specimens, the CCT data and the CT data have approximately 
the same exponent, but the different coefficients mean that 
the CCT requires four times the input of the CT to 
produce identical crack extension. This was not as obvious 
with the aluminum results. 

The regression analysis results for the compact tension 
specimen configurations have been presented graphically in 
figures 3.14, and the fits appear to be good in all cases. 

The aluminum compact tension specimen data were 
investigated further to determine if any thickness effect 
was exhibited. No trend could be seen in the various data, 
and thus it is concluded that no observable thickness effect 
is present. The IN718-STA1 data was treated similarly, but 
it was apparent that B alone defines the nickel CT behavior. 

In summary, the data in tables 3.2 and 3.3 could be 
used for design purposes by matching similar ligament 
configurations . The fit functions were used to determine the 
value of Jjc for each configuration, and the results are 
summarized, below. The Jjc ©stimate is the intersection of 
the power-law fit functions from the tables, above, with a 
line of slope equal to twice the flow stress and offset by 
.008-inches (0.2-millimeters): 

J k = J 0 • ( Aa) m = 2 • c 0 • (Aa-. 008) (3.2) 
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This is the procedure recommended in ASTM E813-87 [3.3]. 
Mathematically, Jj c is derived from an iterative process, 
and the results appear in tables 3.4 and 3.5. 

The values of Jj c varied through the coirplete range of 
ligament configurations of the aluminum compact tension 
specimens, although Rl = 4 and Rl = 2 results were virtually 
identical. This did not agree with the literature [3.4] 
where it had been suggested that Jj c is a constant. The 
regression analysis was generated from relatively extensive 
crack growth, but the results were thought to be bonafide 
and applicable to the determination of Ji c . Too much 
deviation occurs between the extremes to suggest that an 
averaging would be appropriate or that the value is a 
constant. The estimates in this case were derived from J M , 
but the results were virtually the same using the individual 
specimen J-integral results. Applying the same procedures 
to the aluminum CCT specimen results provided different 
approximations, although Jjc f° r the = CCT and the 
side-grooved CCT configurations were the same. 

A look at the Jj c results for the nickel alloy provided 
the same impression (table 3.5). The estimate of Jjc f° r 
the B = 0.80-inch configuration and the side-grooved 
configuration are identical, but the thinner configurations 
have much higher estimates of Jj c - 

The reason for the discrepancy between the data and the 
literature is unknown, but inspection of figures 3.14 
suggests that the differences do exist. It seems reasonable 
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to assume that the extent of the adjustments from the crack 
lengths derived from the calibration functions and the 
unloading compliances to the average physical crack lengths 
may be involved, but the good agreement of the data adjusted 
in this fashion suggests that this would not be the cause of 
the variations of Jjc from configuration to configuration. 

3.3 Conclusio ns and Recommendations 

1. The qualified J^R-curve data of AL6061-T651 compact 
tension specimens followed Rl» the ratio of the 

ligament length-to-ligament thickness. Shorter 
ligament lengths relative to the thickness are 

characterized by lower values of Rl» and these 

configurations provide results that are close to 
plane strain. 

a) A constraint condition of plane strain was 

reached when Rl was reduced to Rj_, ** 1/2. The 
plane strain condition corresponded to the 
lowest and flattest of the J M R-curves and was 
approximated by compact tension specimen tests 
for a side-grooved configuration (10-percent 
per side) . 

b) The condition of plane stress corresponded to 
the highest and steepest of the J M R- curves and 
this was reached by most of the aluminum CCT 
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which appeared to follow the 


specimens, 
blunting line to high values of Jj*. 
c) For aluminum CT specimens, a stable response 
was reached when Rl was increased to Rl *° 2 . 
When Rj_, was increased beyond this, no further 
increase of the Jj^R-curve slope was apparent . 
The J^R- curve generated was lower than that of 
plane stress, and this was associated with the 
gradient of constraint from the surface to the 
center for CT specimens. 

2. For the IN718-STA1, nickel alloy, the basis for 
comparison of the constraint, was the thickness. 

a) The fracture resistance associated with plane 
strain was estimated by testing side-grooved 
compact tension specimens. The constraint 
condition of this greatest thickness (B = 0.85- 
inches) was not equal to plane strain, but was 
slightly higher, indicating a lower constraint 
than plane strain. 

a) The CCT specimens followed the blunting line 
and are taken as the plane stress condition. 

b) The CT specimens with the thinnest ligament 
section provided the steepest J^R-curve, but 
these did not follow the blunting line. 

3. The J M R curve results have been provided in Chapter 
IV, and a regression analysis has been conducted 
with results presented in tables 3.2 and 3.3. The 
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4 . 


power-law functional form produced good agreement 
with the JjvjR-curve data. 

The estimates for Jjq were not constant for the 
variety of configurations tested. ^IC was 

estimated using the ASTM procedure [3.1], and the 
results appear in tables 3.4 and 3.5. 

5. The Jj^R-curve data generated in standard fracture 
tests should be subjected to qualification using a 
validity limits criterion. The best technique was 
the inflection point of the J^R-curve. 

6. JjflR is recommended for the evaluation of fracture 
behavior over the J-integral because correlation is 
extended to a much higher degree of evolution of 
the fracture process. 
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Table 3.1: 




Materials Selected 

AL6061-T651, aluminum alloy 
IN718-STA1, nickel-base super alloy 

Ligament Thickness 

0.25 £ B £ .85 inches 

Ligament Length 

.25 £ b n £ 1.0 inches 

Ligament 

Proportion 

.5 £ b 0 /B £ 4 

Specimen Types 

1/2T and IT compact tension with .5 £a 0 /W£ 

.75, (bending with some tension). 

Center-crack tension, 2W = 2, a n /W = .5, (pure 
tension) , 

Single-edge notch tension, W = 2, a n /W = .75, 

(bending plus tension, intermediate to the other 
two) 

Side-Grooving 

Side-grooved or no 


Conf iquration 


C 


CT. Rr = 2 


C 


SG or R 


14,760 


9379 


3092 


Configuration 


B = 1 /4-inch 


CT. B * 1/2 


C 


C 


SG and B 


10,599 


4232 


2370 


1333 


46,299 


3.2: 
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Figure 3.1: Collections of the JR and J M R behavior of Task I AL6061-T651 and 

IN718-STA1 fracture specimens: (a) JR curve and (b) J M R curve data, all 
aluminum CT specimens, (c) JR curve and (d) J M R curve data, all nickel 
CT specimens. 
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Qualified JR Data for AL6061-T651 Compact Tension Specimens 



• B-1/2 * B-.25, 20% SG 

• B-1/2 20% SG • B-.85 



Crack Extension. Delta-a. Inches 



• B-.25 • B-1/2 • B-.85 • B-.85 


• UB-2 • to®-2 • h/B-4 


Figure 3.2: Collections qualified J resistance curves for AL6061-T65 1 compact tension 

specimens with similar ligament proportions: (a) R L = 1/2 or side-grooved, 
(b) R l = 1, and (c) R L = 2 and R L = 4. 
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Qualified J M R Data for AL6061-T651 Compact Tension Specimens 



. B-1/2 • B-.25. 20% SG 

• B-1/2 20% SG • B-.85 




• B-.25 • B-1/2 • B-.B5 • B-.85 


Figure 3.3: Collections of qualified J M resistance curves for AL6061-T651 compact 

tension specimens with similar ligament proportions: (a) R L = 1/2 or side- 
grooved, (b) R L = 1, and (c) R L = 2 and R L = 4. 
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Qualified J M R Data for AL6061-T651 Compact Tension Specimens 



• b/B=1/2 • b/B=1 • b/B=2,4 • CT.SG 


Figure 3.4: Collection of qualified J M R curves for AL6061-T651 compact tension 

specimens comparing different ligament proportions, R L . 




Qualified J M R Data for IN718-STA1 Compact Tension Specimens 



• t^1/4 • b-1/2 • b-.8 • tvl 



• B-1/2. 6G • B-8, 80 * M 80 



• b-1/4 • b-1/4 • b-1/2 • b-1/2 


Figure 3.5: Collections of qualified J M resistance curves for IN718-STA1 compact 

tension specimens with similar ligament thicknesses or side-grooved: (a) B 
= 1/4, (b) B = 1/2, and (c) side-grooved. 
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J-Modified, lb-in/in A 2 


Qualified J M R Data for IN718-STA1 Compact Tension Specimens 


7000 
6000 
5000 
4000 
3000 
2000 
1000 
0 

0 0.2 0.4 0.6 

Crack Extension, Delta-a, inches 


Figure 3.6: Collection of qualified J M R curves for IN718-STA1 compact tension 

specimens comparing different ligament thicknesses, B. 
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Qualified J M R Data for AL6061-T651 Center-Cracked Tension Specimens 



• B=1/8 • B=1/4 • B=1/2 

• B=.85 • B=.85,SG 


Figure 3.7: Qualified J M R curves for AL6061-T651 center-cracked tension specimens 

comparing different ligament thicknesses, B, (or equivalently ligament 
proportions, RJ with constant ligament length, b = 1/2- inch. 
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Qualified J M R Data for AL6061-T651 Specimens 



Crack Extension, Delta-a. inches 
• B-1/2 • B-.B5 



0 0 2 0.4 

Crack Extension, Deltas, inches 

« SENT.b/B-1/2 • 8ENT4WB-1 


Figure 3.8: Qualified J M R curves for AL6061-T651 specimens: (a) single-edge 

notched tension specimens of different ligament thicknesses, B, (or 
equivalently ligament proportions, R|_) with constant ligament length, b — 
1/2- inch, and (b) comparison of CT specimens with SENT specimens. 
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J-Modified, lb-in/in A 2 


Qualified J M R Data for AL6061-T651 Specimens 


5000 
4000 
3000 
2000 
1000 
0 

0 0.2 0.4 0.6 

Crack Extension, Delta-a, inches 


Figure 3.9: Qualified J M R curves for AL6061-T651 specimens comparing CT, SENT, 

and CCT configurations. 
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J-Modified, lb-in/in A 2 


Qualified J M R Data for IN718-STA1 Center-Cracked Tension Specimens 


7000 
6000 
5000 
4000 
3000 
2000 
1000 
0 

0 0.02 0.04 0.06 

Crack Extension, Delta-a, inches 

• B=1/4 * B=1/4,20%SG 


Figure 3.10: Qualified J M R curves for IN718-STA1 center-cracked tension specimens 
comparing nonside-grooved to side-grooved configurations, B = 1/4-inch. 
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J-Modified, lb-in/in A 2 


Qualified J M R Data for IN718-STA1 Specimens 


7000 
6000 
5000 
4000 
3000 
2000 
1000 
0 

0 0.2 0.4 0.6 

Crack Extension, Delta-a, inches 

Figure 3.11: Q ualifi ed J M R curves for IN718-STA1 specimens comparing CT and CCT 
configurations. 
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Qualified J M R Data for IN718-STA1 Compact Tension Specimens 



• B=.8 • B=1/2, SG • B=.8, SG * B=.8, SG 


Figure 3.12: Collection of qualified J M R curves for IN718-STA1 compact tension 
specimens comparing side-grooved condition to thickest ligament 
configuration, B = 0.8-inches. 
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J-Modified, lb-in/in A 2 


Qualified J M R Data for IN718-STA1 Compact Tension Specimens 

2000 
1500 
1000 
500 
0 

0 0.1 0.2 
Crack Extension, Delta-a, inches 

• B=1/4 • B=1/2 • B=1/2 

Figure 3.13: Collection of qualified J M R curves for IN718-STA1 compact tension 
specimens comparing different ligament proportions, R L . 
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J m R Curves and Regression Results 


AL6061-T651 All CT Specimens 



■ SGARL-1/2 « RL-1 
• RL-2 « RL-4 


IN71 8-STA1 , CT All Thicknesses 



• 66 « B-OJO • B-1/2 « B-1/4 


Figure 3.14: Results of regression analysis, power law fits, (a) AL6061-T65 1 CT 
specimens, and (b) IN718-STA1 CT specimens. 
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